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Abstract We review current ideas on the origin of galactic and extragalactic mag- 
netic fields. We begin by summarizing observations of magnetic fields at cosmological 
redshifts and on cosmological scales. These observations translate into constraints on 
the strength and scale magnetic fields must have during the early stages of galaxy 
formation in order to seed the galactic dynamo. We examine mechanisms for the gen- 
eration of magnetic fields that operate prior during inflation and during subsequent 
phase transitions such as electroweak symmetry breaking and the quark-hadron phase 
transition. The implications of strong primordial magnetic fields for the reionization 
epoch as well as the first generation of stars is discussed in detail. The exotic, early- 
Universe mechanisms are contrasted with astrophysical processes that generate fields 
after recombination. For example, a Biermann-type battery can operate in a proto- 
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galaxy during the early stages of structure formation. Moreover, magnetic fields in 
either an early generation of stars or active galactic nuclei can be dispersed into the 
intergalactic medium. 

Keywords Magnetic fields, Inflation, Early Univers, Quark-Hadron Transition 

There is much to be learned about cosmic magnetic fields. We have a rather 
sketchy information about the field distribution on the largest scales, and the 
origin of the magnetic fields remains a mystery. 

Alexander Vilenkin 2009 



1 Introduction 

Magnetic fields are observed in virtually all astrophysical systems, from planets to 
galaxy clusters. This fact is not surprising since gravitational collapse and gas dynamics, 
the key processes for structure formation, also amplify and maintain magnetic fields. 
Moreover, the conditions necessary for a magnetic dynamo, namely differential rotation 
and turbulence, exist in galaxies which are the building blocks for large scale structure. 
The one notable example where magnetic fields are searched for but not yet found is 
in the surface of last scattering. All this raises an intriguing question: When did the 
first magnetic fields arise? 

Numerous authors have suggested that magnetic fields first appeared in the very 
early Universe. There is strong circumstantial evidence that large scale structure formed 
from the amplification of linear density perturbations that originated as quantum fluc- 
tuations during inflation. It is therefore natural to consider whether quantum fluc- 
tuations in the electromagnetic field might similarly give rise to large-scale magnetic 
fields. 

Indeed, magnetic fields were almost certainly generated during inflation, the elec- 
troweak phase transition, and the quark-hadron phase transition but with what strength 
and on what scale? And more to the point, what happened to these fields as the Uni- 
verse expanded? Were these early Universe fields the seeds for the magnetic fields ob- 
served in present-day galaxies or clusters? And even if not, did they leave an observable 
imprint on the cosmic microwave background? 

Exotic early universe mechanisms for field generation stand in contrast with mech- 
anisms that operate in the post-recombination Universe. There are several ways to 
generate magnetic fields during the epoch of structure formation. At some level, all 
mechanisms begin with a battery, a process that treats positive and negative charges 
differently and thereby drives currents which in turn generate fields. A Biermann bat- 
tery, for example, can (in fact must) operate during the formation of a proto-galaxy. 
While angular momentum in proto-galaxies is generated by the tidal torques due to 
nearby systems, vorticity arises from gasdynamical processes. The same processes also 
drive currents and hence generate fields, albeit of small amplitude. 

The Biermann battery also operates in compact objects such as accretion disks 
and stars. Since the dynamical timescales for these systems is relatively short, tiny 
seed fields are rapidly amplified. The magnetic fields in AGN and/or Pop III stars can 
be expelled into the proto-galactic medium providing another source of seed fields for 
subsequent dynamo action. 
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In this review, we survey ideas on the generation of magnetic fields. Our main focus 
is on mechanisms that operate in the early Universe, either during inflation, or during 
the phase transitions that follow. We contrast these mechanisms with ones that oper- 
ate during the early stages of structure formation though we leave the details of those 
ideas for the subsequent chapter on magnetic fields and the formation of large scale 
structure. The outline of the chapter is as follows: In Section 2, we summarize observa- 
tional evidence for magnetic fields at cosmological redshifts and on supercluster scales 
and beyond. In Section 3, we discuss the generation of magnetic fields during inflation. 
We make the case that inflation is an attractive arena for magnetic-field generation 
but for the fundamental result that electromagnetic fields in the standard Maxwell 
theory and in an expanding, spatially flat, inflating spacetime and massively diluted 
by the expansion of the Universe. However, one can obtain astrophysically interesting 
fields in spatially curved metrics or with non-standard couplings between gravity and 
electromagnetism. Section 4 addresses the question of whether fields can be generated 
during a post-inflation phase transition. We will argue that strong fields almost cer- 
tainly arise but their scales are limited by the Hubble radius at these early times. Only 
through some dynamical process such as an inverse cascade (which requires apprecia- 
ble magnetic helicity) can one obtain astrophysically interesting fields. In Section 5 we 
take, as given, the existence of strong fields from an early Universe phase transitions 
and explore their impact on the post-recombination Universe. In particular, we discuss 
the implications of strong primordial fields on the first generation of stars and on the 
reionization epoch. Finally, in Section 6 we briefly review field-generation mechanisms 
that operate after recombination. A summary and some conclusions are presented in 
Section 7. 



2 Cosmological magnetic fields observed 

The existence of microgauss fields in present-day galaxies and galaxy clusters is well 
established. These fields can be explained by the amplification of small seed fields 
over the 13.7 Gyr history of the Universe. However, there is mounting evidence that 
microgauss fields existed in galaxies when the Universe was a fraction of its present 
age. Moreover, there are hints that magnetic fields exist on supercluster scales. Both of 
these observations present challenges for the seed field hypothesis. In this section, we 
summarize observational evidence for magnetic fields at early times and on cosmological 
scales and briefly discuss the implications for the seed field hypothesis. 



2.1 Galactic magnetic fields at intermediate redshifts 

Microgauss fields are found in present-day galaxies of all types as well as galaxy clusters 

(see, for example. lKronberdl99llWidrowl20ollCarilli and Tavlorl2002l ; lKulsrud and Zweibel 

2008) . Perhaps more significant, at least for our purposes, ar e observations of magnetic 

fields in intermediate and high redshift galaxies. For example. lKronberg. Perry, and Zukowskil 
( 1992) found evidence for a magnetized galaxy at a redshift of z = 0.395. To be spe- 
cific, they mapped the rotation measure (RM) across the absorption-line quasar, PKS 
1229-021 (z = 1.038) and determed the residual rotation measure (RRM - defined to 
be the observed roation measure minus the Galactic rotation measure). The RRM was 
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then identified with an intervening galaxy whose magnetic properties were inferred by 
detailed modelling. 

Along similar lines, IBernet et al. provide a compellin g argument for mag- 

netic fields at z ~ 1.3. Earlier work by Kronberg et al.l (|2008l ~) found a correlation 
between the spread in the quasar RM distribution and their redshift. The naive expec- 
tation is that the spread in the distribution should decrease with redshift. Recall that 
the polarization angle is proportional to the square of the wavelength; the proportion- 
ality constant is what we define as the RM. As electromagnetic radiation propagates 
from source to observer, the rotation angle is preserved by the wavelength increases 
as (1+z) . Hence, RM is diluted by a factor (1 + z) 2 . In principle, the change 
in the RM distribution with red shift could be indica tive of a redshift dependence in 
quasar magnetic fields. However, IBernet et alJ l|2008h sorted the sample according to 
the presence of Mgll absorption lines and showed that the RM spread for set of objects 
with one or more lines was significantly greater than for the objects with no absorption 
lines. Mgll absorption lines arise as the quasar light passes through the halos of nor- 
mal galaxies. The implication is that intervening galaxies produce both large RMs and 
Mgll absorption lines and hence, the intervening galaxies must have strong magnetic 
fields. Simple estimates suggest that the fields are comparable to those in present-day 
gal axies and that these galaxies are at a redshift of z ~ 1.3. 

lAthreva et alJ ( 19981 ) studied 15 radio galaxies with redshifts between z ~ 2 and z — 
3.13 and found significant RM's in almost all of them. Moreover, the RM's were found 
to differ significantly between the two radio lobes suggesting that they are intrinsic to 
the object rather than due to the Faraday screen of the Galaxy. The RM's, corrected 
for cosmological expansion and with the Galactic contribution removed, were typically 
ranged from 100 — 6000 radm~ 2 implying microgauss field strengths. 

Observations of magnetic fields at intermediate redshift imply a shorter time over 
which the dynamo can operate. Consider the standard ylCDM cosmological model 
with n m = 0.272 and fij^ = 0.728 where Q m and Qa are the density, in units 
of the critical density , for matter (both baryons and dark matter) and dark energy 
iKomatsu et~aH (|201of l. In Figured! we show the age of the Universe as a function of 
redshift for this cosmology. We also show the amplication factor for a seed magnetic 
field where we assume exponential growth and one of three choices for the growth rate, 
r = 1.5 Gyr -1 , 2.5 Gyr -1 , or 3.5 Gyr -1 . A seed field of only 10" 21 G is required to 
reach microgauss strength assuming 7~'2.5Gyr _1 . However, for the same growth rate, 
a 10 -11 G seed field is required to reach microguass strengths by a redshift z = 1.3 



2.2 Magnetic fields on supercluster scales and beyond 



iKim et al. I 1119891 ) used the Westerbork Synthesis Radio Telescope to map the Coma 
cluster and its environs at 326 MHz. Their results provide what remains the best direct 
evidence for magnetic fields on supercl uster scales. 

Recently, iNeronov and Vovkl ( 2010h argued that the deficit of GeV gamma rays in 



the direction of TeV gamma-ray sourses yields a lower bound of 3 x 10~ 16 G on the 
strength of intergalactic magnetic fields. The reasoning goes as follows: TeV gamma 
rays and photons from the diffuse extragalactic background light produce pairs 
which, in turn, inverse Compton scatter off photons from the cosmic microwave back- 
ground (CMB). The scattered CMB photons typically have energies in the GeV range. 
In the absence of appreciable magnetic fields, these secondary photons contribute to 
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Fig. 1 Upper pa nel: Age of the Univer se as a function of redshift for the cosmology described 
in the text and in Komatsu ct al. (201CJ). Lower panel: Amplification factor for a seed magnetic 
field assuming exponential growth with one of three choices for the growth rate: F = 1.5 Gyr -1 
r = 2.5 Gyr" 1 or T = 1.5 Gyr" 1 

the overall emission toward the original TeV source. However, magnetic fields will de- 
flect the intermediate pairs. Comparison of model predictions with the observed 
spectrum from HESS Cherenkov Telescopes and upper limits from the NASA F ermi 
Gamma-Ray telescope hint at just such a deficit which [Neronov and Vovkl 1 201Cf ) use 
to derive their lower limit on the magnetic fields. 



3 Magnetic fields from inflation 

3.1 General considerations 

The hierarchical clustering scenario provides a compelling picture for the formation of 
large-scale structure. Linear density perturbations from the early Universe grow via 
gravitational instability. Small-scale objects form first and merge to create systems of 
ever-increasing size. The spectrum of the primordial density perturbations, as inferred 
from the CMB anisotropy spectrum and various statistical measures of large scale 
structure (e.g., the galaxy two-point correlation functi on) is gene r ally t hought to be 
scale- invariant and very close to the form proposed bv IZerdovichl Jl970f). One of the 
great successes of inflation is that it leads to just such a spectrum (|Guth and Pi|[l983 : 
iHawkind Il98l IStarobinskh] 1 19821 V It is therefore natural to ask whether a similar 
mechanism might generate large-scale magnetic fields. 

In order to understand the meaning and significance of the results for density per- 
turbations, we must say a few words about horizons in cosmology. The Hubble radius, 
essentially, the speed of light divided by the Hubble parameter, sets the maximum scale 
over which microphysical processes can operate. In a radiation or matter-dominated 
Universe the Hubble scale is equal, up to constants of order unity, to the causal scale 
which is defined as the distance over which a photon could have propagated since the 
Big Bang. Moreover, the Hubble scale grows linearly with time f in a radiation or 
matter-dominated Universe whereas the physical size of an object associated with a 
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Fig. 2 Evolution of the physical size for the Hubble radius (solid curve) and two scales, 
Ai and A2 (dashed curves) as a function of scale factor a. Shown is the point at which the 
scale Ai crosses outside the Hubble radius during inflation and back inside the Hubble radius 
during the matter-dominated phase. The scale factors ano and a MD correspond to the start 
of the radiation and matter dominated epochs, respectively. Here, Ai enters the Hubble volume 
during the matter-dominated epoch while A2 enters the Hubble volume during the radiation- 
dominated epoch. 

fixed comoving (or present-day) scale grows as t 2 or t 3 / 2 . Therefore, a physical scale 
crosses inside the Hubble radius with smaller objects crossing earlier than larger ones. 
The point is illustrated in Figure [2] and explains why it is so difficult to generate 
large-scale magnetic fields in the early Universe but after inflation. 

During inflation, the Hubble parameter is approximately constant (the spacetime is 
approximately de Sitter) and a physical scale that is initially inside the Hubble radius 
will cross outside the Hubble radius or Hubble scale, at least until some time later in the 
history of the Universe (Figure We therefore have the potential for microphysical 
processes to operate on large scales during inflation with the consequences of these 
processes becoming evident much later when the scale re-enters the Hubble radius. 

Inflation provides the dynamical means for generating density perturbations: quan- 
tum mechanical fluctuations in the de Sitter phase excite modes on the Hubble scale 
with an energy density set by the Hubble parameter. Therefore, to the extent that the 
Hubble parameter is constant during inflation, the energy density in modes as they 
crosses outside the Hubble radius will be scale-independent. 

There is one further and crucial part to the story. During inflation, the energy 
density of the Universe is approximately constant. It is, indeed, the constant energy 
density that drives the exponential expansion of the de Sitter phase. Naively, we expect 
that the energy density in some (relativistic) fluctuation scales as a -4 and is therefore 
diluted by an enormous factor. However, for density perturbations with a coherence 
length greater than the Hubble radius, the energy density scales as a~~ 2 . The ratio, r, 
of the energy density in the fluctuation relative to the background density therefore 
decreases as aT 2 during inflation but grows as a 2 or a during the subsequent radiation 
and matter dominated phases. The net result is that r is the same when the fluctuation 
re-enters the Hubble radius as when it crossed outside the Hubble radius during infla- 
tion. The aT 2 behaviour is often referred to as super-adiabatic growth since it implies 
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that the energy density in the fluctuation grows as a 2 relative to the energy density in 
the radiation field. 

The equation of motion for a scalar field, cf>, in a curved spacetime is 

V Q V Q - - ^ = (1) 

where R is the Ricci scalar, V is the scalar field potential (the driving term for infla- 
tion, but ignored in the present discussion), and £ is a dimensionless constant. Super- 
adiabatic growth occurs for a minimally-coupled scalar field (£ = 0). On the other hand, 
the field evolves adiabatically (energy density scales in the same way as radiation) if it 
is conformally coupled (£ = 1/6). 

The equation of motion for energy density perturbations is identical to that for the 
minimally coupled scalar field. On the other hand, electromagnetism is conformally- 
coupled to gravity (at least in the simplest version of the theory) and any de Sitter- 
induced quantum fluctuations will have an energy density that scales as a~ 4 . One finds 
that the relative strength a magnetic mode at the end of inflation is 




where A is the comoving scale of the mode and rapj ~ 10 GeV is the Planck mass. 
The above ratio also depends on the energy scale of our inflation model (M) and on 
the associated reheating temperature (Tjih). Assuming efficient reheating (i.e. setting 
Trh ~ M), expression ([2]) yields 

^nr 104 A- 4 . (3) 
pt 

During the radiation era, the energy density of the universe is dominated by that 
of the relativistic species (i.e. pt = py oc a~ 4 ) and the high conductivity of the matter 
has been restored. As a result, the magnetic flux is conserved (i.e. B oc a -2 ) and the 
dimensionless ratio 

r= ^l -i° 4 A- 4 , (4) 
Pi 

remains constant. This result implies a field strength no greater than 10~ 50 G on a 
comoving scales of order 10 kpc which are the scales relevant for the galactic dynamo. 

We are lead to the conclusion that inflation-produced magnetic fields are astrophys- 
ically uninteresting. However, this 'negative' result holds for the standard formulation 
of Maxwell's equations and under the assumption of a spatially-flat FLRW cosmology. 
In the next sections, we show that super-adiabatic growth can occur in various "open" 
cosmologies and in models in which certain additional couplings between electromag- 
netic field and gravity are included. 

3.2 Maxwell's equations 

The Maxwell field may be invariantly described by the antisymmetric Faraday tensor, 
F ab . Relative to an observer moving with 4- velocity u a , we can write F a b = 2u[ Q _E b ] + 
e abcB c , where E a = F ab u b and B a = e a t, c F bc /2 represent the electric and magnetic 
components of the EM field respectively. Maxwell's equations then split into two pairs 
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of pr opagation and constraint equations l|Tsagasll200i I Barrow. Maartens. and Tsagasl 
2007). The former consists of 



E(a) =-^&E a + {o ab + uj ab ) E b + e abc A b B c + curlB Q - J a , (5) 



and 



B {a ) = -3 OB a + (a ab + uj ab ) B - e abc A»E c - curlSa , (6) 

which may be seen as the 1+3 covariant analogues of the Ampere and the Faraday 
laws respectively. The constraints, on the other hand, read 

D a E a = fl — 2oj a B a and T) a B a = 2u a E a , (7) 

providing the 1+3 forms of Coulomb's and Gauss' laws respectively. In the above 0, 

a abi ^ab an d A a respectively represent the volume expansi on, the shear, the vorticity 

and t he 4-acceleration associated with the observer's motion ( Tsagas. Challinor. and Maartena 
I2008T ). Also, J a and fi are the 3-current and the charge densities respectively. 

These equations, together with the Einstein equations and the Ricci identities, 
lead to wave equations for the electric and magnetic fields. For example, linearised on 
a Friedmann-Lemaitre-Robertson-Walker (FLRW) background, the wave equation s of 
the electric and the magnetic components of the Maxwell field read (|Tsagasll200"5h 

E a - T) 2 E a = -5HE a + ~ K<J> + 3p)E a - W 2 E a - I TlEa - j a - 3HJ a (8) 

and 

B a - T) 2 B a = -5HB a + I k( P + 3p)B a - 4H 2 B a - ~ TZB a + curlJa , (9) 

respectively. Note that H — a/a = 0/3 is the background Hubble parameter, TZ = 
6K/a 2 represents the Ricci scalar of the spatial sections (with K = 0,±1 being the 
associated 3-curvature index), while n = 8-kG is the gravitational constant. The 3- Ricci 
term arises from the purely geometrical coupling between the electromagnetic field and 
the spacetime curvature. □ We will return to this particular interaction to examine the 
way it can affect the evolution of cosmological magnetic fields. 



3.3 Ohm's law in the expanding Universe 

The literature contains various expressions of Ohm's law, which provides the prop- 
agation equation of the electric 3-current. For a single fluid at the limit of resistive 
magnetohydro dynamics (MHD), Ohm's law takes the simple form I Greenberd Il97l1 : 
|jacksonlll975l ) 

Ja = sE a , (10) 

with ? representing the electric conductivity of the matter. In highly conducting en- 
vironments, c — > 00 and the electric field vanishes. This is the familiar ideal-MHD 
approximation where the electric currents keep the magnetic field frozen-in with the 
charged fluid. Conversely, when the conductivity is very low, q — > 0. Then, the 3- 
currents vanish despite the presence of nonzero electric fields. Here, we will consider 
these two limiting cases. For any intermediate case, one needs a model for the electrical 
conductivity of the cosmic medium. 



1 For the f ull expression s of Eqs. J5) and 10, written in a general spacetime, the reader 
is referred to [Tsagas (200^). There, one can also see how the different parts of the geometry 
(i.e. the Ricci and the Weyl fields) can affect the propagation of electromagnetic signals. 
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3.4 Adiabatic decay of magnetic fields in a spatially flat FLRW cosmology 

Consider the case of a poorly conductive environment s where there are no 3-currents. 
The wave equation, @, then reduces to (Tsagas 2005) 

B a - B 2 B a = -5HB a + i(p + 3p)B a - 4H 2 B a - 1l ab B b . (11) 

To simplify the above, we introduce the rescaled the magnetic field B a = a, 2 B a and 
employ conformal time, rj, where r\ = 1/a. Then, on using the harmonic splitting 
B a = B( n )Qa n ^ ~ so that D £?(„) = 0, expression (QTJ takes the compact form 

B' ( ' n) +n 2 B (n) = -2KB (n) , (12) 

with the primes denoting conformal-time derivatives and K = 0, ±1 (|Tsagasl[2005l ). 
Note the magneto-curvature term in the right-hand side of (|12|l , which shows that the 
magnetic evolution also depends on the spatial geometry of the FLRW spacetime. 

When the background has Euclidean spatial hyp ersurf aces, the 3-curvature index 
is zero (i.e. K = 0) and expression (|12[) assumes the Minkowski-like form 

B'( n) + n 2 B (n) = . (13) 

This equation accepts the oscillatory solution B( n \ = C\ sin(rw;) + Ci sm(nr\) , which 
recasts into 

B (n) = [Ci sin(nj7) + C 2 cos^)] , (14) 

for the actual B-field. In other words, the adiabatic (Br n ) oc a~ 2 ) depletion of the 
magnetic component is guaranteed, provided the background spacetime is a spatially 
flat and the electrical conductivity remains very poor. This result leads directly to 
Equation [21 

The adiabatic decay-law also holds in highly conductive environments. There, c — > 
oo and, according to Ohm's law (see Eq. (|10|) the electric field vanishes in the frame of 
the fluid. As a result, Faraday's law (see Eq. (J6]l) linearises to 

B a = —IBB a , (15) 

around an FLRW background. The above ensures that B a oc a~ 2 on all scales, regard- 
less of the equation of state of the matter and of the background 3-curvature. This 
result leads directly to Equation [4] 



3.5 Superadiabatic magnetic amplification in spatially open FLRW models 

The "negative" results discussed at the beginning of this section have been largely 
attributed to the conformal invariance of Maxwell's equations and to the conformal 
flatness of the Friedmannian spacetimes. The two are thought to guarantee an adiabatic 
decay-rate for all large-scale magnetic fields at all times. Solution (|14[) . however, only 
holds for the spatially flat FLRW cosmology. Although all three FLRW universes are 
conformally flat, only the spatially flat model is globally conformal to Minkowski space. 
For the rest, the conformal mappings are local. Put another way, in spatially curved 
Friedmann universes, the conformal factor is no longer the cosmological scale factor 
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but has an additional spatial dependence l|Stefanilll990l : lKeane and Barrettl lioOO). The 
wave equation of the rescaled magnetic field (B a = a 2 B a ) takes the simple Minkowski- 
like form (|13p only on FLRW backgrounds with zero 3-curvature. In any other case, 
there is an additional curvature-related term (see expressions and Ijl2| l1. which 
reflects the non-Euclidean spatial geometry of the host spacetime. As a result, when 
linearised around an FLRW background with nonzero spatial curvature, t he magnetic 
wave equation reads ( Tsagas and K andus 20051; iBarrow and Tsagadliool ) 



£(„) + 



(n 2 ± 2) 



#(„) = , 



(16) 



with the plus and the minus signs indicating the spatially closed and the spatially open 
model respectively. Recall that in the former case the eigenvalues are discrete (with 
n 2 > 3), while in the latter continuous (with n 2 > 0). As expected, in either case, the 
curvature-related effects fade away as we move on to successively smaller scales (i.e. for 
n 2 » 2). 

Following (|16[) . on FLRW backgrounds with spherical spatial hypersurfaces, the 
B-field still decays adiabatically. The picture changes when the background FLRW 
model is open. There, the hyperbolic geometry of the 3-D hypersurfaces alters the 
nature of the magnetic wave equation on large enough scales (i.e. when < n 2 < 2). 
These wavelengths include what we may regard as the largest subcurvature modes 
(i.e. those with 1 < n 2 < 2) and the supercurvature lengths (having < n 2 < 1). 
Note that eigenvalues w ith n 2 = 1 correspond to the curvature scale with physical 
wavelength A = Xx — a ( Lvth and WoszczvnalE995l ) . Here, we will focus on the largest 
subcurvature m odes. 

In line with iTsaeas and Kandus! l|2005l ) and IBarrow and Tsagasl fcoosh . we intro- 



duce the scale-parameter k 2 = 2 — n 2 , with < k 2 < 2. Then, fc z = 1 indicates the 
curvature scale, the range < k 2 < 1 corresponds to the largest subcurvature modes 



and their supercurvature counterparts are contained with in the 1 < k 
In the new notation and with K = — 1, Eq. (|16|) reads 



< 2 interval. 



0. 



(17) 



with the solution given by &(k) = C\ sinh(|fc|?7) + C% cosh(|fe|rj). Written in terms of 
the actual magnetic field, the latter takes the form 



B 



(fc) 



Cie 



|fc| (»7— ^o) 



+ C 2 e 



(-) 

Van/ 



(18) 



Magnetic fields that obey the above can experience superadiabatic amplification with- 
out modifying conventional electromagnetism and despite the conformal flatness of the 
FLRW host. For instance, during the radiation era, the scale factor of an open FLRW 
universe evolves as a oc sinh(n). Focusing on the curvature length, for simplicity, we 
may set |fc| = 1 in Eq. (|18[) . On that scale, the dominant magnetic mode never decays 
faster than Br % \ oc a" 1 . The B-field has been superadiabatically amplified[3 



2 Mathematically, the most straightforward case of superadiabatic amplification occurs on 
a Milne background. The latter corresponds to an empty spacetime with hyperbolic spatial 
geometry and can be used to describe a low density open universe. The scale factor of the 
Milne universe evolves as a oc e" 7 , which substituted into solution H18H leads to 



B(k) — ^5 



(?) 



\k\-2 



-C e 



-1*1-2 



(19) 
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Analogous amplification also occurs in open FLRW univers es with an inflationary 

(i.e. p = — p) equation of state. There, the scale factor evolves as I Tsagas. Challinor. and MaartenU 



2008) 



where 17,770 < 0. Substituting the above into Eq. (|18|l . we find that near the curvature 
scale (i.e. for |fc| — > 1) the magnetic evolution is given by 

with C3, C4 depending on the initial conditions. This result also implies a superadiabatic 
type of amplification for the B-field, since the dominant magnetic mode never decays 
faster than B/jj oc a -1 . The adiabatic decay rate is only recovered at the end of 
inflation, as r\ — > 00 

The strength of the residual magnetic field is calculated in a way analogous to that 
given in the previous section. In particular, when B oc a~\ we find that 



by the end of inflation ( Tsagas and Kandusl 1200*51 ; iBarrow and Tsagaslliooi ). As in 



Eq. (|2|, M and Trh are measured in GeV, while A is given in Mpc. Thus, the higher the 
scale of inflation, the lower the reheatint temperature and also the smaller the number 
of e-folds, the stronger the amplification. Setting M ~ 10 17 GeV and Trh ~ 10 9 GeV, 
for example, the current (comoving) magnetic strength varies between ~ 10 -35 and 
~ 10~ 33 Gauss on scales close to the present curvature scale. The latter lies between 
~ 10 4 and ~ 10 5 Mpc when 1 — Q ~ 10~ 2 today. We note that the condition 1 — 
Q ~ 10~ 2 is within t he values allowed by recent analysis of the WMAP observations 
(jKomatsu et al.l [2010h These lengths are far larger than 10 kpc; the minimum magnetic 
size required for the dynamo to work. Nevertheless, once the galaxy formation starts, 
the field lines should break up and reconnect on scales similar to that of the collapsing 
protogalactic cloud. Note that the above quoted strengths assume that the ratio r — 
PbIPi remains constant after inflation. As we have seen earlier, however, magnetic 
fields spanning lengths close to the curvature scale are superadiabatically amplified 
during the radiation era as well. When this is also taken into account, the residual 
B-field increases further and can reach magnitudes of up to 10~ 16 G at present. 

Galactic-scale magnetic fields of strength 10 -16 are stronger than those generated 
by many of the other scenarios considered in the literature (see below) and may be 
strong enough to seed the galactic dynamo. Recall, however, that inflation was intro- 
duced to avoid various shortcomings of the standard cosmological model including the 



Consequently, all magnetic modes spanning scales with < k 2 < 2 are superadiabatically 
amplified. Close to the curvature scale, that is for k 2 — > 1, the dominant magnetic mode is 
-B(!) oc a -1 . Stronger amplification is achieved on supercurvature lengths, with Bffe) oc a^ 2 ^ 2 
at the homogeneous limit (i.e. as k 2 2). 

3 The magnetogeometrical interaction and the resulting effects are possible because, when 
applied to spatially curved FLRW models, inflation does not lead to a globally flat de Sitter 
space. Although the inflationary expansion dramatically increases the curvature radius of the 
universe, it does not change its spatial geometry. Unless the universe was perfectly flat from 
the beginning, there is always a scale where the 3-curvature effects are important. It is on 
these lengths that primordial B-fields can be superadiabatically amplified. 
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so-called flatness problem. Essentially, inflation is meant to inflate away (push to ex- 
tremely large scales) any curvature that might exist in the pre-inflationary Universe. 
Evidently, for superadiabatic growth of magnetic field without explicit conformal sym- 
metry breaking (as described in the next section) one requires enough inflation to push 
the curvature scale beyond our present Hubble volume, but n ot too far beyond. To 
be quantitative, the most recent analysis by the WMAP group (jKomatsu et alJl2O10h 
finds that the effective energy density parameter for curvature, = 1 — Sl^ — i? m , 
is constrained to be —0.0133 < < 0.0084. Should future measurements find J?j. 
to be inconsistent with zero and positive, they would lend credence to the idea of 
superadiabatic magnetic amplification by the mechanism described above. 



3.6 Inflation-produced magnetic fields via non-conformal couplings 

As pointed out above, electromagnetism is conformally coupled to gravity and there- 
fore, in a spatially flat, FLRW cosmology, the magnetic fields generated during in- 
flation decay adiaba t ically and are therefore of negligible astrophysical importance. 
iTurner and Widrowl (|l988T ) pointed out that by adding additional terms to the La- 
grangian such as RA 2 and R^vA^A" one explicitly breaks conformal invariance and 
can essentially force the mag netic field to behave l i ke a m inimally-coupled scalar field. 

The terms introduced bv lTurner and Widrowl (|l988T ) also break gauge mvariance 
and hence induce an effective mass for the photon whose size depends on the spacetime 
curvature. In fact, the effective mass-squared is negative for the case where the magnetic 
field behaves as a minimally-coupled scalar field. The negative mass-squared signals an 
instability in the semi-classical equations for the field and can be viewed as the origin 
of the superadiabaticity. The se terms also give rise to ghosts in the theory w hich signal 



in instability of the vacuum jHimmetoglu. Contaldi. and Pelosd l2"o09 2009). A theory 



with ghosts is internally consistent on ly as an effective low-energy theory and h ence is 
only valid below some energy scale, A. Illimmetoglu. Contaldi. and Pelosol 1 20091 ) argue 
that A MeV but this scale is well below the scales assumed in models for inflation- 
produced magnetic fields. Thus, their results call into question the viability of the 
mechanism. 

Numerous authors have attempte d to find more effect i ve an d more natural ways 
to break conformal invariance. Indeed. [Turner and Widrowl 1 1988h . recognizing the po- 
tential difficulties associated with the iL4 2 -terms con sidered terms such as RF 2 which 
break conformal invariance but not gauge invariance. iRatral |l992) demonstrated that 
appreciable magnetic fields couple be produced during inflation if the electromagnetic 
field couples to the inflaton field, $ , through a term of the form e a ^ ' F illV F il ' v where a is a 
constant. In his model, the inflaton potential is also an exponential: V (4>) oc exp (— 
An attractive feature of this model is that is preserves gauge invariance since the ad- 
ditional term is const ructed from the Maxwell tensor, F, ra t her t han the gauge field, 
A. Along similar lines iGasperini. Giovannini. and Venezianol 1 19951 ) demonstrated that 
magnetic fields of sufficient strength to seed the galactic dynamo could be produced 
in a string-inspired cosmology. In their model, the electromagnetic field is coupled to 
the dilaton which, in turn, is coupled to gravity. The dilaton is a scalar field that nat- 
urally arises in theories with extra dimensions and whose vacuum expectation value 
effectively controls Newton's constant. A detailed and critical a nalysis of magnetic-field 
produ ction in string-inspired m odels of inflation ca n be found in lMartin and Yokovamal 
(2008) (see, also, the review bv lSubram anian 201o]). There, particular attention is paid 
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Fig. 3 The cosmological periods. Starting with reheating the quantum chromodynamic quark- 
gluon-plasma period lasts until hadron formation. Non-Abelian Weibel instabilities may be 
responsible for generating seed magnetic fields during this time which in presence of free 
energy grow from non-Abelian thermal fluctuations. At later times the period of relativistic and 
classical plasma sets on , Here seed magnetic fields can be generated from thermal background 
electrodynamic fluctutations if thermal anisotropies or beams can exist in this regime. 

to a potential back-reaction problem which arises if the electric fields produced along 
with the magnetic fields have an energy density comparable to the background energy 
density. 

4 Magnetic fields from early universe phase transitions 

The early universe was characterized by a series of phase transitions in which the nature 
of particles and fields changed in fundamental ways (see, Figure [3]). For example, elec- 
troweak symmetry breaking marked the transition from a high-energy regime in which 
the W and Z bosons and the photon were effectively massless and interchangeable to 
one in which the W and Z bosons were heavy while the photon remained massless. The 
transition also marked the emergence of two distinct forces: electromagnetism and the 
weak nuclear force. Likewise, the Quark-Hadron phase transition marked the transition 
from the free quark-gluon phase to one in which quarks were locked into baryons. Both 
of these transitions have the potential to generate strong magnetic fields since they in- 
volve the release of an enormous amount of free energy and since they involve charged 
particles which can, in turn, drive currents. Indeed, strong magnetic fields are almost 
certainly generated. The question is one of physical scale since the Hubble radius is so 
small at these early times. 

4.1 General considerations 

At any phase in the history of the Universe, the strength of the fields generated by some 
microphysical process is limited by equipartition with the background energy density. 
Moreover, the maximum scale for magnetic fields generated by a microphysical process 
is set by the Hubble radius. For a radiation-dominated Universe, the energy density 
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of the Universe is given by p oc #*T 4 where T is the temperature of the thermal bath 
and g* counts the effective number o f relativistic degrees of freedom (see, for example, 
iKolb and Turnerlll990l ; |Peacocklll999h . The Hubble radius is given by L H =c/H where 
H = a~ 1 da/dt oc T 2 is the Hubble parameter. Numerically, we have 

Bmax (I = H~\ T) = SequipartitioB CO K 10 18 Gauss ( 15Q ^ eV ) ' (23) 



for 



^lOOcm(l^) 2 . (24) 

In practive, the strength of fields generated during some early Universe phase transition 
will be well below the value set by equipartition and have a length scale significantly 
smaller than the Hubble radius. We therefore write B(l, T) = fB max for I = gc/H and 
/ and g are constants. 

We are typically interested in fields o n scales much larger than the Hubble scale 
of some early Universe phase transition. iHoganl |l983h argued on purely geometric 



grounds, that the field strength on a large scale L due to small-scale cells of size I with 
field strength B(l) will be B(L) = B(l) (l/L) 3 ^ 2 . Moreover, in the absence of some 
dynamical effect, the field strength will be diluted by the expansion as a -2 . Thus, a 
field B[l, T) generated prior to recombination will lead to a field at recombination with 
strength 

B(L, T roc ) = B(l, T) (^2£) 2 * . (25) 

More generally we have the following scaling: 

B (L, T) ~ f g 3 / 2 T~ 3/2 L~ 3/2 . (26) 



4.2 First-order phase transitions 

Detailed calculations of magnetic field generation during the electroweak and QCD 
phase transtions have been carried out by numerous authors. By and large, these 
groups assume that the transitions are first-order, that is characterized by a mixed- 
phase regime in which bubbles of the new phase nucleate and expand, eventually filling 
the v olume. The energy associated wi t h the bubble walls is released as a form of latent 
heat. IQuashnock. Loeb. and Spergell (1989) demonstrated that a Biermann battery 
can operate during the QCD phase transition. The up, down, and strange quarks (the 
three lightest quarks) have charges 2/3, -1/3, and -1/3 respectively. If these quarks 
were equal in mass, the quark-gluon plasma would be electrically neutral. However, 
the strange quark is heavier and therefore less abundant. The implication is that there 
is a net positive charge in the quark-gluon plasma and a net negative charge in the 
lepton sector. Electric currents are therefore generated a t the bubble walls that separate 
the qu ark phase from the baryon phase sweep through. IQuashnock. Loeb. and Spergell 
(|l989l ) found that 5 G fields could be generated on scales of 100 cm at the time of the 
QCD phase transition. Following the arguments outlined above, the field strength on 

galactic scales at the time of recombination wou ld be (a disappointly sma ll) ~ 10~ 31 G. 

So mewhat larger estimates were obtained bv lCheng and Olintol l|l994h and lSigl. Olinto. and Jedamzikl 
(1997) who realized that as the hadronic regio ns grow, baryons would concentra te on 



the bubble walls due to a "snowplow" effect. IjSigl. Olinto. and Jedamzikl ( |l997h also 
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showed that fluid instabilities could give rise to strong magnetic fields during the QCD 
phase transition.) For reasonable param eters, they obtained fields about seven orders 
of magnitude larger than those found bv lQuashnock, Loeb. and Spergel 1 19891 ). 

Magnetic fields can arise during cosmological phase transitions even if they are sec- 
ond order — that is, ph ase transitio n s sign aled by the smooth and continuous transition 
of an order parameter. IVachaspatil jl99ll ). for example, showed that gradients in the 
Higgs field vacuum expectation value (the order parameter for the electroweak phase 
transition) induce magnetic fields on a scale ~ ^ew w ^ n strength of order Q^^V^EW 
where T^w ls the temperature of the electroweak phase transitio n and Qrw is the 
Higgs field coupling constant. To estimate the field on larger scales, IVachaspatil (|l99lh 
assumed that the Higgs field expectation value executed a random walk with step size 
equal to the original coherence length. The field strengths were small (10~ 23 G on 
100 kpc scales) but not neglibible. 



4.3 Inverse cascade 

The discussion above suggests that strong magnetic fields are likely to have been gen- 
erated in the early Universe but that their coherence length is so small, the effective 
large-scale fields are inconsequential for astrophysics. However, dynamical mechanisms 
may lead to an increase in the coherence length of magnetic fields produced at early 
times. Chief among these is an inverse cascade of magnetic energy from small to larg e 
scales which occurs when there is substantial magnetic helicity l|Frisch et al.l Il975l ) . 
The e ffect was investig at ed in the c o ntext of primordia l magn e tic fields (see, for ex - 
ample. ICornwalll dl997h : [Sor] (|l999T ): iField and Carrolll l|2000h : iBrandenburj (jl996l ): 



iBaneriee and Jedamzikl l|2004h . As the Universe expands, magnetic energy shifts to 
large scales as the field attempts to achieve equilib r ium w hile conserving magnetic he- 
licity. Under suitable conditions. IField and Carrolll (2000) showed that astrophysically 
interesting fields with strength 10 -10 G could be generated on 10 kpc scales. 



4.4 Plasma processes capable of generating magnetic fields 
4-4- 1 Chromodynamic magnetic fields? 

The QCD regime lasts from the end of reheating until hadron formation, roughly 
tns ~ 10~ 6 s after the Big Bang (see Figure [3j). In this regime, matter comprises 
massive bosons, gluons, quarks and leptons and forms a hot dense chromoplasma or 
quark-giuon plasma (QGP). At the higher temperatures, that is, not too long after 
reheating, the QGP is asymptotically free and can be considered collisionless. Since 
many of the particles in the QGP carry electric charge, under certain conditions, they 
can generate induced Yang-Mills currents ja(x) — D^F^ V (x), with Yang- Mills field 
Fftv = A u ,^ — A/x,!/ — ig[A/x, A v \ expressed through the non-Abelian gauge field A a ;v{x). 
The colour index a corresponds to the N 2 — 1 colour channels. These currents couple 
to the electromagnetic gauge field and consequently produce magnetic fields. Several 
mechanisms associated with the QCD phase transition were discussed in the previous 
section. Here we explore whether thermal plasma instabilities in the QGP can lead to 
appreciable fields. 
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The simplest plasma mecha nism is the W eibel (current filamentationjf] instability 
discovered in classical plasma (Weibel 1959). Its free energy is provided by a local 
pressure anisotropy A — Pn/P±_ — 1 7^ in the non-magnetic plasma. ||,_L refer to 

the two orthogonal directions ||,_L of t he pressure tensor P = P±\ + (Pn — Pj_)|||| 
whose non-diagonal elements are small ( Blaizot and Ianucll2002h . Pressure anisotropy 



creates microscopic currents and hence microscopic magnetic fields. Free energy can 
also be provided by partonic beams passing the QGP. Such beams naturally introduce 
a preferred direction may cause additional pressure anisotropy by dissipating their 
momentum in some (collisionless) way. 



merical simulations ( Arnold and Moore 2006a 


2006b; 


Arnold 2007a; Arnold and Moore! 


2007b|; Arnold and Leang|2007c 


; Rebhan. Strickland, and Attemsj2008l; Romatschke and Rebhanl 


200d: Ischenkell2008l: Istrickland 


l200rJ. 20071) Drove that the OCD beam-driven Weibel 



instability excites magnetic fields which subsequently scatter and thermalize the beams 
by magnetising the partons. The linear waves, that is, oscillations of the effective quark 
phase space momentum distribution $ eJ y(p) as function of the 4- momentum, p, are so- 
lutions of the sem i -class ical QGP dispersion relation i Pokrovskv and Selikhovi ri988: 
lMr6wczvrisklll98ci Il993h in Fourier space k — (oj/c, k) 



detpc 2 ^ - kk? 



2 in 
- lo e 



(|k|)] = 



(27) 



with permeability (velocity v l = p 1 / \J p l p\ is the velocity) given as functional of the 
effective phase space density $ eJ y(p), which does not depend any more on the colour 
index 



'(w,k) 



5 ij , g_ 



d 3 P v l {d4> eff (p)/dp 

8tt 3 w - k ■ v + iO 



ii r 



(uj — k • v)<5 + k l v 



(28) 



Instability is found at low frequency, uj ~ 0, non-oscillatory filamentation modes with 
wave vectors, kj_, perpendicular to the beam four- velocity, U. The implication is that 
stationary magnetic fields are generated. 

Analytical growth rates of transverse modes, where Imcj > 0, have been obtained 
for simple gaussian and ot h er moc k equilibrium particle distributions and nuclear 
physics parameters. lArnoldl ( 2007ah has shown that the breakdown of perturbation 
theory at momenta p ~ g 2 T and the fact that theory becomes non-perturbative at 
high T implies that it can be treated as if one had T — plus weak coupling. Thus, 
the semi-classical approach describes long-range properties In numerical simulations 
one takes adva ntage of this fact, linear ises around a stationary homogeneous locally 
colourless state ( Blaizot and Ianucll2002h . and considers the evolution of the fluctuation 
W^(v, x) of the distribution function according to the non-Abelian Vlasov equation and 
Yang-Mills current density: 



v„D fi W v (v,x) = -VaF av (x), f{x) = -g- 



d 3 p p" d$ e ff(p) 



8tt 3 



dp a 



(29) 



4 For a discussion of its physics in classical plasmas sec Fried (1959). 

5 The scalar potential A picks up a Debye screening mass mu and decouples at distances ^> 
m D ~ (fl-^ 1 ) - 1 leaving the vector potential (transverse chromo-electromagnetic) fields, which 
is equivalent to classical plasmas where at frequencies uj > u) v above the plasma frequency uj p 
any propagating perturbation is purely electromagnetic. 
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Fig. 4 Numeric al simulation results of We ibel fields under Abelian and non-Abelian (QGP) 
conditions (after Arnold and Moore 2006a). The Abelian and spatially ld-cases have a long 
linear phase of exponential growth. The corresponding linear phase in the spatially 3d-case is 
very short (shaded region), followed by a nonlinearly growing phase which reaches maximum 
field strength and afterward increases weakly and linearly. Closer investigation has shown that 
this further linear growth is caused by cascading to shorter wavelength (larger k) related to 
increasing frequencies thereby dispersing the free energy that feeds the nonlinear phase into 
a broad spectrum of magnetic oscillations. The blue horizontal line is the nonlinear satura- 
tion level of the low-frequency long-wavelength Weibel magnetic field that persists during the 
cascade and presumably survives it. 

where v v = (l,p). The latter enters the field equations which close the system and 
describes the evolution of fields and particles, 



- 9 (E + vxB)-V p * eJ j 
D„ /•■''"' = -g 2 (2Tr)- 3 f^pv" (DT^E + v x B) ■ V p $ e 



(30) 



'/i r ' — ~ y V^" > J u yu yj-j j ynj t v /v jj; ■ v p-^ eJ Ff 

Instab ility requires that the effective distribution $ eJ y (p) is anisotropic ( Arnold and Moorel 
2006a) in p to lowest order 0(1). The maximum unstable wave vector k m and growth 
rate 7 m = +(Im 

iS)rn of the QCD Weibel mode scale as 



In 



*eff(p)/|p| 



(31) 



p> 



where moo ~ 9\/N>/p> is the "effective mass" scale defined by the spatial number 
density N> of the particles with momentum p> that contribute to the integral and 
anisotropy, i.e. N> = $ eJ y (p)d 3 p/87r 3 . 

Figure 3] plots three simulation results: an Abelian run, a spatially Id, and a spa- 
tially 3d-non- Abelian run. Shown is the magnetic fluctuation energy density B 2 /8n as 
function of time (all in proper simulation units). The Abelian and the ld-non- Abelian 
cases cover just their linear (exponentially growing) phases. The 3d-non- Abelian case 
differs from these in several important respects. Its linear phase is very short, shown as 
the shaded region. It is followed by a longer nonlinear growth phase when the magnetic 
energy density increases at a slower and time dependent rate until reaching maximum, 
when it starts decaying. Afterwards it recovers to end up in a further slow but now 
purely algebraic linear growth. 
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Weibel saturation level. The nonlinear phase results from the nonlinear (wave-particle 
interaction) terms that come progressively into play when the field energy increases. 
The subsequent brief decay phase and the following slow linear growth phase result from 
the sudden onset of a (turbulent) long-wavelength cascade to shor ter wavelength and 



higher frequencies. This result h as been convincingly demonstrated ( Arnold and Moora 
l2006bl : lArnold and Leanell2007ch . The cascade transfers the excess energy that the non- 
linear instability feeds per time unit into the long wavelengths to shorter wavelengths 
in such a way that in the average the energy density in the long wavelengths (low 
frequencies) stays at a constant level which, presumably, survives the cascade. This 
value, for the settings of the simulation, is 

t-Weibel ~ 3m »/j 2 (32) 



For the coupling constant one may use l|Bethke!l2009h the "world-average value" g 



Ana s (M z ) « An x 0.12 = 0.48vr, where M z is the mass of the Z boson, while for 
additional knowledge is required of the effective parton distribution $ eJ y , i.e. the state 
of the undisturbed distribution including the anisotropy. 

Unfortunately, the above formula cannot be easily used to estimate the long- 
wavelength Weibel field saturation value in the early universe. If we accept that the 
saturation level is robust within few orders of magnitude, then, because m 2 ^ ~ ojp /V0, 
one has as for an estimate (in physical units) 



B sat ttLu% v / 2fi h/c 3 6 (33) 

where u) p is the effective chromo-plasma frequency uip = gT/mj^Xjj which is related to 
the temperature T, Debye mass mp, and screening length Xjj, and 6 ~ tan -1 (A -1 ) 
is an effective anisotropy angle ( Arnold and Leans! l2007cn . With these numbers one 



estimates quite a strong saturation magnetic field 

B sat » KT 19 (m e /M z Ve)n eff ~2x W~ 25 n eff /V8 G (34) 

a value that depends linearly on the effective parton density n e ff. Taking, say, n e g ~ 
10 10 cm -3 , it yields a large QCD saturated low- frequency seed magnetic field of the 
order of B sat ~ 10 _E VG. This seems unrealistically high. What it shows, however, 
is that only very small anisotropies 9 ~0(1) are required for production of magnetic 
fields during the QGP phase even if the above estimate is wrong by several orders of 
magnitude. It seems thus worthy of further considering this possibility of primordial 
magnetic field generation. 

The mass density in the standard model varies by about 40 orders of magni- 
tude from reheating until quark-hadron transition where it is njf' « 10 20 kg/m 3 . 
At weak unification it is of the order of n™ w 10 33 kg/m 3 . Number densities de- 
pend on the dominant particle mass chosen. If, for simplicity, we assume a mass 
M ~ 10 3 GeV, then number densities at the quark-hadron transition are still of the 
order of n^J ~ 10 21 cm -3 . With such high effective densities one obtains extreme 

values for the saturated mag netic field the order of B sat ~ 10" 4 G and even larger at 
an earlier time. This is probably unrealistic and thus could hardly be right. Clearly, 
if such high fields existed they must have been attenuated during further evolution 
(expansion) of the universe to the low values needed after recombination ended. 
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Fig. 5 The spectrum S(k) of magnetic fields in the pre- recombination era generated from 
cosmological density perturbations (after Ichiki et al. 2006) plotted in units of magnetic field 
(G) and shown to be composed of contributions from baryon-phon slip and anisotropic photon 
stresses. Below roughly k <1/Mpc the baryon slip dominates the spectrum while for larger k 
(shorter wavelength than roughly 10 Mpc) the anisotropic stresses contribute most. 



Thermal fluctuation level. The argument might be weakened when no anisotropy 
exists at all. Then one is led to the determination of the thermal magnetic Weibel fluc- 
tuation level which, of course, is well below the saturation level estimated above. For 
its estimation the effective isotropic equilibrium distribution $ eJ y(p) and the complex 
response function, i.e. the QGP polarisation tensor U iiV , are needed whose determi- 
nation requires solution of the Vlasov equation. A rough estimate of the thermal level 
can be found from the above simulation results by taking advantage of the evolution 
equation of the average magnetic energy density (B 2 (t)) in long wavelength fluctua- 
tions and assuming that the magnetic energy density measured at a certain time t\ in 
the linear phase only evolved from the thermal level {b 2 (t = 0)} at time t = according 
to 

<B 2 (*i)) ~ (b 2 (t = 0)} exp 2 7m t (35) 

^From the linear phase in Figure [4] one finds that 7 m as 0.28/moo- This value used 
in the last equation yields an approximate (initial) thermal Weibel magnetic energy 
density level of (b 2 ) w 2.88 x W^^m^g^ 2 , which corresponds to a thermal Weibel 
magnetic fluctuation of average amplitude 

(bweibel) ~ 3-4 x HT 29 n e # G (36) 



during the corresponding QCD phase; at the quark-hadron transition this becomes 
(bf ht ) ~100 fiG, which still is very large for thermal background fields. Levels such high 
as the estimate would, probably, be subject to attenuation in the classical plasma phase 
before decoupling if they should account for the current large scale fields. Of course, 
this value can only give a hint on the possible importance of thermal fluctuations as it 
has been taken and rescaled from the available simulations which have been performed 
for other purposes (thermalisation of nuclear matter beams). Whether realistic or not 
can be decided only after developing a QCD theory for the cosmological QGP phase. 



20 



4- 4-2 Thermal fluctuations shortly before photon-matter decoupling. 

The previous section dealt with the possible generation of low frequency seed magnetic 
fields during the QCD ph ase of the early univers e. A more conventional proposal has 
been elaborated recently ( Ichiki et al.| [2006, 2007) and is proposed to work during the 



classical plasma phase before recombination. These authors make the reasonable as- 
sumption that during this epoch the coupling between photons and electrons by Comp- 
ton scattering is much stronger than the coupling between photons and ions. Cautious 
examination of the particle-photon interaction in a cosmological density fluctuation 
field indeed shows that pressure anisotropy and currents are induced. The generalised 
Ohm's law that includes the photon interaction allows for a finite electron current flow 
because of the differences in the bulk electron and proton velocities caused. This effect 
produces the desired magnetic fields. However, only the second order density perturba- 
tions in the Compton scattering terms lead to fields that survive on a range of spatial 
scales. The power spectrum S(k) of the fields in the long wavelength range (see Figure 
[5]) scales as \J k 3 S(k) oc k Here the photon-caused anisotropic stress dominates. 

Seed magnetic fields reach values of Bsiml0~ ls G on scales of ~1 Mpc and B ~ 
10~ 14 G on ~10 kpc scales. After decoupling these fields decay adiabatically with 
expansion of the universe. In a standard cosmology they should today be of strength 
B(to) ~ 10 -24 G at 1 Mpc and ~ 10~ 20 G at 10 kpc and may have played a role in 
structure formation after recombination. 



5 Magnetic fields and the cosmic microwave background 

The field of cosmology has witnessed a revolution lead, in large part, by detailed mea- 
surements of the CMB anisotropy and polarization spectra. Fundamental cosmological 
parameters such density parameters of baryons, dark matter, and dark energy and the 
Hubble constant are now known to a precision unimaginable just two decades ago. 

If magnetic fields were present at the time of matter-radiation decoupling or soon 
after, then they wou ld have an effect on the anisotropy and polarization of the CMB (see 
ISubramanianl (|2006h : TPurrerl l|2007h for reviews). First, a very large scale (effectively 
homogeneous) field would select out a special direction, lead to anisotropic expansion 
around this direction, hence leading t o a quadrupole anisotropy in the cosmic microwave 
background (CMB) (see, for example. iThornell 1967h . The degree of isotropy of the CMB 
then i mplies a limit of several nG on the strength of such a field redshifted to the current 
epoch (jBarrow. Ferreira and Suklll997h . 

Primordial magnetogenesis scenarios on the other hand generally lead to tangled 
fields, plausibly Gaussian random, characterized by say a spectrum M(k). This spec- 
trum is normalized by giving the field strength Bq, at some fiducial scale, and as mea- 
sured at the present epoch, assuming it decreases with expansion as B = Bo/a (t), 
where a(t) is the expansion factor. Since magnetic and radiation energy densities both 
scale with expansion as 1/a 4 , we can characterize the magnetic field effect by the ratio 
Bo/(8ttp 7 o) ~ 10~ 7 B 2 9 where p 7 is the present day energy density in radiation, and 
_B_9 = Bq/(10~ s G). Magnetic stresses are therefore small compared to the radiation 
pressure for nano Gauss fields. 

Nevertheless, the scalar, vector and tensor parts of the perturbed stress tensor 
associated with primordial magnetic fields lead to corresponding metric perturbations, 
including gravitational waves. Further the compressible part of the Lorentz force leads 
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to compressible (scalar) fluid velocity and associated density perturbations, while its 
vortical part leads to vortical (vector) fluid velocity perturbation. These magnetically 
induced metric and velocity perturbations lead to both large and small angular scale 
anisotropics in the CMB temperature and polarization. 

The scalar contribution has been th e most subtle to calculate, and has only be- 
gun to be understood by several groups dGiovannini and Kunzj [ 2008l ; lYamazaki et al.l 
I2OO8I : iFinelli. Paci and Paolettlll2008l : IShaw and Lewij201(J l." The anisotropic stress as- 
sociated with the magnetic field leads in particular to the possibility of two types of 
scalar modes, a potentially dominant mode which arises before neutrino decoupling 
sourced by the magnetic anisotropic stress. And a compensated mode which remains 
after the growing neutrino a nsiotr oic stress has compensated the magnetic anisotropic 
stress (cf. Shaw and Lewis! (l2010h for detailed discussion). The magnetically induced 
compressible fluid perturbations, also changes to the acoustic peak struc ture of the an- 
gular anisotropy power spectrum (see, for example. I Adams et al.l [l996). However, for 
nano Gauss fields, the CMB anisotropics due to the magnetized scalar mode are grossly 
subdominant to the anisotropics generated by scalar perturbations of the inflaton. 

Potentially more important is the co ntribution of the Alfven mode driven by the ro- 
tational component of the Lorentz force (S ubramanian a nd BarrowH T998; Mac k. Kahniashvili and Kosowskvl 
|2002| ; ISubramanian. Seshadri and BarrowlbOOsTnLewislboO^) . Unlike the compressional 
mode, whi ch gets strongly damped below the Silk scale, Lg due to radiative viscosity 
(Silk 1968), the Alfven mode behaves like an over damped oscillator. This is basically 
because the phase velocity of oscillations, in this case the Alfven velocity, is Va ~ 
3.8 x 10 _4 C-B_g much smaller than the relativistic sound speed c/y/3. Note that for an 
over damped oscillator there is one normal mode which is strongly damped and another 
where the velocity starts from zero and freezes at the terminal velocity till the damp- 
ing becomes weak at a latter epoch. The net result is that the Alfven mode survives 
Silk damping d own to much smaller scales; La ~ (Va/c)L c ; <C L,s, the canonical Silk 
damp ing scale Uedamzik. Katalinic and Olinto et al.lll998l : ISubramanian and Barrowl 
1998). The resulting baryon velocity leads to a CMB temperature anisotropy, AT ~ 
5^,K(B-g/3) 2 for a scale invariant spectrum, peaked below the Silk damping scale 
(angular wavenumbers I > 10 3 ). 

The magnetic anisotropic stress also induces tensor perturbations, resulting in a 
comparable CM B temperature anisotropy, but now peak ed on large angular scales of a 
degree or more ( Durrer. Ferreira and Kahniashvilil 120071 ). Both the vector and tensor 

perturbations lead to ten tim es smaller B-type polarization anisotropy, at respectively 

small and large angular scales jSeshadri and Subramanian|200ll ; ISubramanian. Seshadri and Barrowl 
l2003l : lMack. Kahniashvili and Kosowskvll2002l ; lLewij2004l l. Note that inflationary gen- 
erated scalar perturbations only produce the E-type mode. The small angular scale vec- 
tor cont ribution in particular can potentially help to isolate the magnetically induced 
signals I Subramanian. Seshadri and Barrowll2003l ). 

A crucial difference between the magnetically induced CMB anisotropy signals 
compared to those induced by inflationary scalar and tensor perturbations, concerns 
the statistics associated with the signals. Primordial magnetic fields lead to non- 
Gaussian statistics of the CMB anisotropics even at the lowest order, as magnetic 
stresses and the temperature anisotropy they induce depend quadratically on the 
magnetic field. In contrast, CMB non-Gaussianity due to inflationary scalar pertur- 
bations arises only as a higher order eff ect. A computation of the nongaussianity of the 
magnetically induced signal ha s begun jSeshadri and Subramanianll2009l ; ICaprini et all 
l2009l ; ICai. Hu and Zhan5l2010l) . based on earlier calculations of non-Gaussianity in the 
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magnetic stress energy ( Brown and Crittenden! [20051 ). This new direction of research 
promises to lead to tighter constraints or a detection of strong enough primordial mag- 
netic fields. 

A primordial magnetic field leads to a number of other effects on the CMB which 
can probe its existence: Such a field in the inter galactic medium can cause Faraday 
rotation of the polarized component of the CMB, l eading to the generation o f new 
B-type signals from the inflationary E-mode signal i Kosowskv and Loebl 1199(1 ). Any 
large-scale helical component of the field leads to a parity violation eff ect, inducing 
non-zero T-B and E-B cross-correlations ( Kahniashvili and Ratrall2005l) : such cross- 
correlations between signals of even and odd parity are neccessarily zero in standard 
inflationary models. The damping of primordial fields in the pre-reombination era can 
lead to spectral distortions of the CMB (Jedamzik, Katalinic & Olinto 2000), while their 
damping in the post-recombination era can change the ionization and thermal history 
of t he universe and hence the electr o n scattering optical depth as a f u nction of redshift 
(see Sethi and Subramanianl ( 20051 ); iTashiro and Sugivamal 1 2006al ); ISchleicher et al.l 
( 20081 ). and Section 5.3). Future CMB probes like PLANCK can potentially detect the 
modified CMB anisotropy signal from such partial re-ionization. In summary primordial 
magnetic fields of a few nG lead to a rich variety of effects on the CMB and thus are 
potentially detectable via observation of CMB anisotropies. 



6 Implications of strong primordial fields in the post-recombination 
universe 

If strong magnetic fields have been produced during phase-transitions in the early 
universe, and if these fields had some non-zero helicity, they may have remained 



strong until recombination and be yond ( Christensson. Hindmarsh. and Brandenburd 
l200ll : iBaneriee and Jedamzikll20oi ). They could then affect the thermal and chemical 
evolution during the dark ages of the Universe, the formation of the first stars, and the 
epoch of reionization. 



6.1 Implications during the dark ages 

At high redshift z > 40, the universe is close to homogeneous, and the evolution of 
the temperature, T, is governed by the competition of adiabatic cooling, Compton 
scattering with the CMB and, in the presence of strong magnetic fields, ambipolar 
diffusion. It is thus given as 

dT _ 8o- T a fl r r 4 ad 3: e (T - T rad ) 
dz 3H(z)(l + z)m e c 1 + fn e + x e 
2T _ 2(£ AD - L coo i) 
1 + 2 3nk B H{z)(l + z)' 1 ' 

where Lad i s the heating functio n due to ambipolar diffusion (AD), L coo \ the cool- 
ing function ( Anninos et al.l fl997t ) . o"t the Thomson scattering cross section, cir the 



Stefan-Boltzmann radiation constant, m e the electron mass, c the speed of light, kg 
Boltzmann's constant, n the total number density, x e = n c /njj the electron fraction 
per hydrogen atom, T rax ± the CMB temperature, H(z) is the Hubble factor and /jj e is 
the number ratio of He and H nuclei. 
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AD occurs due to the friction between ionized and neutral species, as only the 
former are directly coupled to the magnetic field. Primordial gas consists of several 
neutral and ionized species, and f or an appropriate description of th is process, we thus 
adopt the multi-fluid approach of IPinto. Galli. and Bacciottil 1 20081 '). defining the AD 
heating rate as 

i AD = ! |^l(VxB)xB/ J B| 2 , (38) 

where ?7ad is given as 

^AD=E^AD,„- (39) 

n 

In this expression, the sum includes all neutral species n, and ??AD,n denotes the AD 
resistivity of the neutral species n. We note that the AD resistivities themselves are a 
function of magnetic field strength, temperature and chemical composition. 

In the primordial IGM, the dominant contributions to the total resistivity are 
the resistivities of atomic hydrogen and helium due to collisi ons with protons . Thes e 
are calculated based on the momentum transfer coefficients of IPinto and Gaul (|2008h . 
As the power spectrum of the magnetic field is unknown, we estimate the expres- 
sion in Eq. (|38|) base d on the coherence length Lr , given as th e characteristic scale 
for Alfven damping ([je damzik, Katalinic and Olinto et al.l 1 19981 : ISubramanianl Il998l : 
ISeshadri and SubramanianlfeoOlT ). Contributions from decaying MHD turbulence may 
also b e considered, but are negligible compared to the AD heating ( Sethi and Subramanianl 
l2005h . 

The additional heat input provided by AD affects the evolution of the ionized 
fraction of hydrogen, x p , which is given as 

dx p = [x e x p n H a H - P H {1 - x p )e~ h v VH ^l kT ] 
dz H(z)(l + z)[l + K H (A H + p H )n H (l - x p )] 

x [1 + K H A H n H (l x p )] jggf^J - (40) 

Here, njj is the number density of hydrogen at oms and ions, hp Planck's constant, 
fei on is the collisional ionization rate coefficient IjAbel et al ] |l997l ). Further details of 
notation, as well as the parametrized case B recombinatio n coefficient for atomic hy- 
drogen aji, are given by |S ea ger. Sasselov. and Scottl (| 19991 ). The chemical evolution of 
the primordial gas is solved with a system of rate equations for the chemical species 
HP, H^", H2, HeH + , D, D + , D~, HD + an d HD based on the primordial rate coeffi- 
cients tabulated by ISchleicher et al (2008). For the mutual neutralization rate of H 
and H + , we use the more recent result of IStenrup. Larson, and Elanderl |2009:). The 
evolution of the magnetic energy density Eg = B 2 /8tv is given as 

dE B 4 dp Eg 

— = 3mT~ AD - (41) 

The first term describes the evolution of the magnetic field in a homogeneous universe 
in the absence of specific magnetic energy generation or dissipation mechanisms. The 
second term accounts for corrections due to energy dissipation via AD. 

The dynamical implications of magnetic fields can be assessed from the magnetic 
Jeans mass, the critical mass scale for gravitati onal forces to overcome magnetic pres- 
sure. In the large-scale IGM, it is given as ( Subramanian|[l998l ; ISethi and Subramanianl 
120051 ) 



Mj ~ 10 10 M Q 
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Fig. 6 The evolution of the comoving magnetic field strength due to AD as a function of 
redshift for different initial comoving field strengths, from the homogeneous medium at z = 
1300 to virialization at z = 20. 



Due to ambipolar diffusion, strong magnetic fields also affect the gas temperature and 
thus the thermal Jeans mass, i.e. the critical mass scale required to overcome gas 
pressure. It is defined as 

MTrpN- 1 / 2 / 5k B T \ 3/2 



Mj = hrJ ) (43) 

with Boltzmann's constant kg, the mean molecular weight p and the proton mass m p . 

To d escribe v i rializa tion in the first minihalos, we employ the spherical collapse 
model of lPeeblesI (| 19931 ) for pressureless dark matter until an overdensity of ~ 200 is 
reached. Equating cosmic time with the timescale from the spherical collapse model 
allows one to calculate the overdensity p/pb as a function of time or redshift. In this 
model, we further assume that the formation of the protocloud will reduce the coherence 
length of the magnetic field to the size of the cloud. 

The evolution of the magnetic field strength, the IGM te mperature and the chem - 



ical abundances of different species have been calculated by Schleicher et al. (2009b) 



using an extension of the recombination code RECFAST I Seager. Sasselov. and ScottI 
Il999l ). The results are shown in Figs. 1 6181 As shown in Fig. (0, ambipolar diffusion 
primarily affects magnetic fields with initial comoving field strengths of 0.2 nG or less. 
For stronger fields, the dissipation of only a small fraction of their energy increases 
the temperature and the ionization fraction of the IGM to such an extent that AD 
becomes less effective. For comoving field strengths up to ~ 0.1 nG, the additional 
heat from ambipolar diffusion is rather modest and the gas in the IGM cools below the 
CMB temperature due to adiabatic expansion. However, it can increase significantly 
for stronger fields and reaches ~ 10 4 K for a comoving field strength of 1 nG, where 
Lyman a cooling and collisional ionization become efficient and prevent a further in- 
crease in temperature. The increased temperature enhances the ionization fraction and 
leads to larger molecule abundances at the onset of star formation. 
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Fig. 7 The gas temperature evolution in the IGM as a function of redshift for different 
comoving field strengths, from the homogeneous medium at z = 1300 to virialization at z = 20. 
For the case with Bq = 0.01 nG, we find no difference in the thermal evolution compared to 
the zero-field case. 
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Fig. 8 The evolution of ionization degree, H2 and HD abundances as a function of redshift for 
different comoving field strengths, from the homogeneous medium at 2 = 1300 to virialization 
at z = 20. For the case with Bo = 0.01 nG, we find no difference in the chemical evolution 
compared to the zero-field case. 



6.2 Implications for the formation of the first stars 



We now explore in more detail the consequences of magnetic fields for the formation 
of the first stars, during the protostellar collapse phase. For this purpose, a model 
de scribing the chemical and thermal evolutio n during free-fall c ollapse was developed 
by Idover and Savi n| (|2009h and extended by ISchleicher et~aTl (|2009bh for the effects 
of magnetic fields. Particularly important with respect to this application is the fact 
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that it correctly models the evolution of the ionization degree and the transition at 
densities of ~ 10 8 cm -3 where Li -1 " beco mes the main charge carrier. Based on the 
Larson- Penston type self-similar solution ( Larson! 1 19691 ; |Penstonlll969l ; lYahil 1983), we 



evaluate how the collapse timescale is affected by the thermodynamics of the gas. 

During protostellar collapse, magnetic fields are typically found to scale as a power- 
law with density p. Assuming ideal MHD with flux freezing and spherical symmetry, 
one expects a scaling with p 2 / 3 in the case of weak fields. Deviations from spherical 
symmetry such as e xpected for dynamically imp ortant fields give rise to shallower seal 



ings, e.g. B oc o ' 6 IIBaneriee and Pudritj 12006). B oc p 1 / 2 1 Hennebelle and Fromang 



Machida et al 



l2008l:)Hennebelle and TeyssieifeoolT Based on numerical simulations of 
(|2006r i. we find an empirical scale law 

( M \ °- 0116 

« = 0.57(^J ■ (44) 

In a collapsing cloud, the more general expression for the magnetic Jeans mass, 

Mj = — (45) 

is adopted. In this prescription, $ = ivr 2 B denotes the magnetic flux, G the gravita- 
tional constant, r an appropriate length scale. The calculation of the magnetic Jeans 
mass thus requires an assumption regarding the size of the dense region. Numeri- 
cal hydrodyn amics simulations show that they are usua lly comparable to the thermal 
Jeans length ( Abel et al.ll2002 l; lBromm anc ~Larson||2004l ). This is also suggested by an- 



alytical models for gravitational collapse (|Larsonlll969l : lPenstonlll9 69: Yahil 1983). To 



account for magnetic energy dissipation via AD, we calculate the AD heating rate from 
Eq. (|38|) and correct the magnetic field strength accordingly. We note that due to the 
large range of densities during protostellar collapse, additional processes need to be 
taken into account to calculate the AD resistivity correctly. In particular, at a density 
of ~ 10 9 cm~ 3 , the three-body H2 formation rates start to increase the H2 abun- 
dance significantly, such that the gas is fully molecular at densities of ~ 10 11 cm~ 3 . 
As a further complication, the proton abundance drops con siderably at densities of 
10 8 cm" 3 , such that Li + becomes the main charge carrier ( Maki and Haiimell2004l : 



iGlover and Savinll2009l ). These effects are incorporated in our multi-fluid approach. 

As an initial condition for these model calculations, we use the physical field 
strength and the chemical abundances obtained from the spherical collapse model in 
section lrTTl The relation between co-moving and physical field strength at the beginning 
of this calculation is thus given in Table [1] 



B [nG] B [nG] 

1 4.8 x 10 a 

0.3 6.3 x 10 2 

0.1 1.3 x 10 2 

0.01 1.0 x 10 1 



Table 1 The physical field strength B at beginning of collapse as a function of the comoving 
field strength Bq used to initialize the IGM calculation at z = 1300. 



27 



10000 



1000 



100 




10" W 10* 10" 10" 
Number density [cm~ 



10 1 



10 1 



Fig. 9 The gas temperature as a function of density for different comoving field strengths. 
For Bq = 0.01 nG, the thermal evolution corresponds to the zero-field case. 
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Fig. 10 Thermal (thin lines) and magnetic (thick lines) Jeans mass as a function of density 
for different comoving field strengths. For Bo = 0.01 nG, the thermal evolution and thus the 
thermal Jeans mass corresponds to the zero-field case. 



Fig.[5]shows the calculated temperature evolution as a function of density for differ- 
ent comoving field strengths. For comoving fields of 0.01 nG or less, there is virtually no 
difference in the temperature evolution from the zero-field case. For comoving fields of 
~ 0.1 nG, cooling wins over the additional heat input in the early phase of collapse, and 
the temperature decreases slightly below the zero-field value at densities of 10 3 cm~ 3 . 
At higher densities, the additional heat input dominates over cooling and the temper- 
ature steadily increases. At densities of ~ 10 9 cm -3 , the abundance of protons drops 
considerably and increases the AD resistivity defined in Eq. (|39[1 and the heating rate 
until Li + becomes the main charge carrier. In particular for comoving fields larger than 
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~ 0.1 nG, this transition is reflected by a small bump in the temperature evolution 
due to the increased heating rate in this density range. 

Apart from the transition where Li + becomes the dominant charge carrier, the 
magnetic field strength usually increases more rapidly than p ' 5 , and weak fields in- 
crease more rapidly than strong fields. This is what one naively expects from Eq. (|44l) . 
and it is not significantly affected by magnetic energy dissipation. Another important 
point is that comoving fields of only 10" 5 nG are amplified to values of ~ 1 nG at a 
density of 10 3 cm~ 3 . Such fields are req uired to drive protostellar outflows that can 
magnetize the IGM ijMachida et alJl2006h . 

Fig. llOl shows the evolution of the thermal and magnetic Jeans mass during collapse. 
The thermal Jeans masses are quite different initially, but as the temperatures reach 
the same order of magnitude during collapse, the same holds for the thermal Jeans 
mass. The thermal Jeans mass in this late phase has only a weak dependence on the 
field strength. As expected, the magnetic Jeans masses are much more sensitive to the 
magnetic field strength, and initially differ by about two orders of magnitude for one 
order of magnitude difference in the field strength. For comoving fields of ~ 1 nG, the 
magnetic Jeans mass dominates over the thermal one and thus determines the mass 
scale of the protocloud. For ~ 0.3 nG, both masses are roughly comparable, while 
for weaker fields the thermal Jeans mass dominates. The magnetic Jeans mass shows 
features both due to magnetic energy dissipation, but also due to a change in the 
thermal Jeans mass, which sets the typical length scale and thus the magnetic flux in 
the case that Mj > Mj . 

The uncertainties in these models have been explored further by ISchleicher et"al] 
(|2009ph . an d an independent calculation inclu ding stronger magnetic fields has been 
provided bv lSethi. Haiman. and Pandevll|2010h . We also note that the consequences of 
initially weak magnetic fields for primordial star formation are explored in more detail 
in the next chapter. 



6.3 Implications for reionization 

Strong magnetic fields may influence the epoch of reionization in various ways. As dis- 
cussed above, they may affect the formation of the first stars and change their mass 
scale, and thus their feedback effects concerning cosmic reionization and metal enrich- 
ment. They may further give rise to fluctuations in t he large-scale density field and 
potentially enhance high-redshift structure formation ([Kim. Olinto. and Rosner| [l996: 
ISethi and Subramanianl2005l : lTashiro and Sugivamall2006ah . On the other hand, the in- 
creased thermal and magnetic pressure may indeed suppress star formation in small ha- 
los an d thus delay reionization ( Schleicher et al]|2008l ; lRodrigues. de Souza. and Opherl 



l2010l ). In both cases, unique signatures of the magnetic field may become inprinted in 
the 21 cm signature of reionization, whic h may help to constrain or detect such magnetic 
fields with LOFAE0, EDGEsQor SKA0 (|Tashiro and Sugivamal2006bl ; ISchleicher. Baneriee. and Klessenl 



l2009al ). Upcoming observations of these facilities may thus provide a unique opportu- 
nity to probe high-redshift magnetic fields in more detail. 



6 LOFAR homepage: 

7 EDGES homepage: 



http://www.lofar.org/! 



http://www.haystack.mit.edu/ast/arrays/Edges/ 
SKA homepage: |http: / /www.skatelescope.org/| 
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7 Seed fields in the post-recombination Universe 

We have seen that magnetic fields arise naturally during inflation and during phase 
transitions after inflation but before recombination. However, the effective field strength 
on galactic scales may well be exceedingly small. Indeed, the seed fields for the galactic 
dynamo may well arise from astrophysical processes rather than exotic early-Universe 
ones. In this section, we review three processes which can generate magnetic field in 
the post-recombination Universe. 



7.1 Biermann battery 

In the hierarchical clustering scenario, proto-galaxies acquire angular momentum from 
tidal torques produced by neighboring systems ( Hovldll949l ; |Peebleslll969l : I Whitelll984l ) . 



However, these purely gravitational forces cannot generate vorticity (the gravitational 
force can be written as the gradient of a potential whose curl is identially zero) and 
therefore the existence of vorticity must arise from gasdynamical processes such as 
those that occur in shocks. More specifically, vorticity is generated whenever one has 
crossed pressure and density gradients. In an ionized plasma, this situation drives cur- 
rents which, in turn, generate magnetic field. This mecha nism, known as t he Biermann 
battery and originally studied in th e context of stars JBierman nl 1950T) was consid - 



ered in the cosmological c onte xt by Pudritz and Silkl (|1989T ): Kulsrud et al.l fll997l ') 
Davies and Widrowl (feooch and |x u et al.l (|2008T ). A simple order of magnitude esti 



mate yields 

^biermann ^ ~ 3 X 1CT 21 ^ ^-1^-1 ^ Gauss (46) 

where u> is the vorticity. Since the vorticity in the solar neighborhood is of order 
30kms _1 kpc -1 we expect seed fields of order 10~ 19 G. 



7.2 First-generation stars 

The first generation of stars provides another potential source of seed fields for the 
galactic dynamo. Even if stars are born without magnetic fields, a Biermann battery 
will generate weak fields which can then be rapidly amplified by a stellar dynamo. If 
the star subsequently explodes or loses a significant amount of mass through stellar 
winds, the fields will find their way into the interste l lar me dium and spread throughout 
the (proto) galaxy. Simple estimates bv lSvrovatskiil 1 1970h illustrate the viability of the 
mechanism. There have been some 10 8 supernovae over the lifetime of the galaxy, each 
of which spreads material through a (10 pc) 3 volume. Using values for the field strength 
typical of the Crab nebula, one therefore expects the galaxy to be filled by 10 pc regions 
with field strengths ~ 3 (iG. Assuming the same L -3 / 2 scaling discussed above, one 
finds a field strength of 10~ n G on lOkpc scales, a value significantly larger than the 
ones obtained by more exotic early Universe mechanisms. 

The strong amplification of seed magnetic fields during primordi al star formation 
has b e en suggested in a number of w orks. Analytical estima tes by IPudritz and Silkl 
( 19891 ): iTan and Blackmanl $004) and lSilk and Langeij l|2006l ) suggest that large-scale 



dynamos as well as the magneto-rotational instability could significantly amplify weak 
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magnetic seed fields until saturation. Even before, during the protostellar collapse 
phase, the small-scale dynamo leads to an exponenti al growth of the magnetic fields, 
as fou nd in both semi-analytic and numerical studies (|Schleicher et aljboict ISur et all 

boid ). 



7.3 Active galactic nuclei 

Strong magnetic fields almost certainly arise in active galactic nucleii (AGN). These 
fields will find their way into the intergalactic mediu m via jets thereby providing a 
potential source of m agnetic field for normal galaxies l|Hovlelll969l : lRees!ll987l . 1 19941 : 



iDalv and Loeblll990h . The potential field strengths due to this mechanism may be 
estimated as follows: The rotational energy associated with the central compact (mass 
M) object which powers the AGN can be parametrized as fMc 2 where / < 1. If we 
assume equipartition between rotational and magnetic energy within a central volume 
V c , we find 



If this field then expands adiabatically to fill a "galactic" volume V g one finds B g = 
B c (y./yJ 2/3 . lHovlel(ll969h considered values M = 10 9 M , / = 0.1, V g ~ (100 kpc) 3 
and found B c ~ 10 9 G and B g ~ 10~ 5 G. 



8 Conclusions 

The origin of the seed magnetic fields necessary to prime the galactic dynamo remains 
a mystery and one that has become more, rather than less, perplexing over the years as 
observations have pushed back the epoch of microgauss galactic fields to a time when 
the Universe was a third its present age. Numerous authors have explored the possible 
that the galactic fields observed today and at intermediate redshift have their origin 
in the very early Universe. 

The impetus for the study of early Universe magnetic fields came from the suc- 
cessful marriage of ideas from particle physics and cosmology that occurred during the 
latter half of the last century. It is a remarkable prediction of modern cosmology that 
the particles and fields of the present-day Universe emerged during phase transitions 
a fraction of a second after the Big Bang. The electromagnetic and weak interactions 
became distinct during the electroweak phase transitions at t ~ 10~ 12 s while baryons 
replaced the quark-gluon plasma at t ~ 10~ 6 s. Perhaps more fantastical is the notion 
that galaxies, clusters, and superclusters arose from quantum-produced density per- 
turbations that originated during inflation at even earlier times. This idea is supported 
by strong circumstantial from the CMB anisotropy spectrum, so much so, that it is 
now part of the standard lore of modern cosmology. 

Both inflation and early Universe phase transitions have many of the ingredients 
necessary for the creation of magnetic fields. If our understanding of inflation-produced 
density perturbations is correct, the similar quantum fluctuations in the electromag- 
netic field will lead to fields on the scales of galaxies and beyond. As well, electro- 
magnetic currents, and hence fields, will almost certainly be driven during both the 
electroweak and quark-hadron phase transitions. 



31 



Early Universe schemes for magnetic field generation, however, face serious chal- 
lenges. In the standard electromagnetic theory, and in a flat or closed FLRW cosmology, 
inflation-produced fields at diluted by the expansion to utterly negligible levels. One 
may address this issue by considering fields in an open Universe and "just-so" inflation 
scenario, that is, a scenario with just enough e-folds of inflation to solve the flatness 
problem. Alternatively, one may incorporate additional couplings of the field to gravity 
into the theory though many terms lead to unwanted consequences which render the 
theory unphysical. Observations and advances in theoretical physics may settle the is- 
sue. For example, if future determinations find that the Universe has a slight negative 
curvature (density parameter for matter and dark energy slightly less than one), it 
would give some credance to the idea of superadiabatic field amplification in an open 
Universe. On the other hand, string theory may point us to couplings between gravity 
and electromagnetism that naturally generate fields during inflation. 

The main difficulty with fields generated from phase transitions arises from the 
small Hubble scale in the very early Universe. Strong fields are almost certainly pro- 
duced by one of a number of mechanism. But the coherence length is so small, the 
effective field strength on galactic scales is likely to be well-below the level of interest 
for astrophysics. An inverse cascade may help; if the field has a net helicity, the mag- 
netic field energy will be efficiently transferred from small to large scales. But even if 
fields are uninteresting for the galactic dynamo, they may have an effect on processes 
in the post-recombination Universe such as reionization and the formation of the first 
generation of stars. 

Where, if not the early Universe, did the seed fields for the galactic dynamo arise? 
Astrophysics provides a number of promising alternatives. Galactic disks have angular 
momentum and vorticity. While the former is generated by the gravitational interaction 
between neighboring protogalaxies, the latter comes form gasdynamical effects which 
necessarily generate magnetic fields via the Biermann battery. As well, magnetic fields, 
rapidly created and amplified inside some early generation of stars or in active galactic 
nucleii, can be dispersed throughout the intergalactic medium. 
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